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Abstract. 

We study the conformational properties of polymers in presence of extended 
columnar defects of parallel orientation. Two classes of macromolecules are considered: 
the so-called partially directed polymers with preferred orientation along direction 
of the external stretching field and semiflexible polymers. We are working within 
the frames of lattice models: partially directed self-avoiding walks (PDSAWs) and 
biased self-avoiding walks (BSAWs). Our numerical analysis of PDSAWs reveals, that 
competition between the stretching field and anisotropy caused by presence of extended 
defects leads to existing of three characteristic length scales in the system. At each fixed 
concentration of disorder we found a transition point, where the influence of extended 
defects is exactly counterbalanced by the stretching field. Numerical simulations 
of BSAWs in anisotropic environment reveal an increase of polymer stiffness. In 
particular, the persistence length of semiflexible polymers increases in presence of 
disorder. 
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1. Introduction 

Conformational properties of polymers are the subject of great interest during decades 
due to their important role both in nature and industry. In statistical description of 
macromolecules it is established, that they are characterized by a number of properties, 
which are universal, i.e. independent on any details of microscopic chemical structure 
HIE]. In particular, the mean-squared end-to-end distance (i?^) of long flexible polymer 
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chain in a good solvent scales with number of monomers (molecular weight) N according 
to: 


{R^) ~ 


( 1 ) 


here (...) denotes averaging over an ensemble of conformations, and the universal size 
exponent z/ = 0.587597(7) in three dimensions [3]. 

Many synthetic polymers with a carbon backbone, such as polyethylene, can be 
considered as long flexible polymer chains: the molecular weight of repeating units 
(monomers) is rather small and neighboring carbon atoms can easily rotate around single 
bonds. However, not all polymers satisfy this criterion. In particular, the important 
organic macromolecules such as DNA and some proteins contain massive repeating 
units (e.g., radicals of amino acids) with complicated interactions between them, which 
thus can not freely rotate on their bonds. This class of polymers is usually referred 
to as semiflexible m- An important characteristic here is a persistence length Ip-. 
when the total chain length is much larger than Ip, any polymer behaves as a flexible 
chain obeying the scaling law ([I]), whereas for the chain length much smaller than the 
persistence length, the polymer attains the limit of rigid rod with (i?^) ~ A^T This type 
of “coil-to-rod” conformational transition can be describe by presenting a size measure 
of polymer in a scaling form [6]: 

~ Np (0 . ( 2 ) 

Here, f{x) is a scaling function: 


const, N « Ip, 

x2u-2 N» Ip. 


( 3 ) 


Usually, the semiflexible polymers are characterized by a relatively large value of 
persistence length Ip (e.g., the mechanical persistence length of a DNA macromolecule is 
of order of 50 nm |5]). Conformational properties of semiflexible polymers are intensively 
studied [MS]. 

Another interesting class of polymer macromolecules are so-called directed polymers 
mi. Presence of anisotropy in the system, e.g. the external electrical held in the 
solution of polyelectrolytes, causes the stretching of polymer chain in direction parallel 
to the held, but allows huctuations in transverse directions. Due to anisotropy of 
directed polymers, one should distinguish between the transverse {R\) and longitudinal 
(7?y) components of the end-to-end distance, which scale with the molecular weight N 
according to: 


(i?5) ~ AT^^ii, {Rl)r^N^''^ (4) 

where v\\ = 1 and n_L = 1/2 in all dimensions [T8] . 

Universal scaling properties of polymers can be successfully described within the 
framework of lattice model of self-avoiding random walk (SAW) [19]. The fact, that 
the trajectory can not cross itself rehects the excluded volume ehect in polymers. This 
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model can be generalized for the case of semiflexible polymers by introducing different 
statistical weights for “trans” step (taken in the same direction as the previous one) and 
“gauche” step (leading to a bending of the trajectory), which is known in literature as 
biased self-avoiding walks (BSAW) p[^[20] . The persistence length of BSAW trajectory 
is thus introduced as inversely proportional to the probability of a gauche step. The 
properties of directed polymers can be described within a model of directed random walk 
(DSAW). The DSAW trajectory is considered to be directed along some coordinate axis 
(say x), if projection of each step of the trajectory onto this axis is non-negative [TTjlTQ]. 
One can treat this as a result of stretching held acting in direction x. Note that the 
self-avoidance of directed walk trajectory (excluded volume effect of polymers) does not 
play any role, and the size exponents in (jl]) are trivial in all dimensions [T8] . 

In real physical systems, one often encounters various forms of structural 
inhomogeneities and irregularities. In studying the polymer systems, of great importance 
is understanding of the behavior of macromolecules in the presence of structural 
obstacles, e.g. in colloidal solutions EZH, in gels or intra- and extracellular 
environment [2^12^ . One can introduce disorder into the lattice model of polymer 
by considering some sites of the lattice as occupied by defects and forbidden for SAW 
trajectory. 

In general, dealing with disordered systems one usually encounters two types 
of ensemble averaging, treated as annealed and quenched disorder [25ll26]. In the 
former case, a characteristic time of impurities dynamics is comparable to relaxation 
times in the pure system, whereas in the latter case the impurities are hxed and the 
conhgurational average over an ensemble of disordered systems with different realization 
of the disorder is needed. Though in general the critical behavior of systems with 
quenched and annealed disorder could be quite different, an equivalence of the influence 
of quenched and annealed disorder on scaling properties of long flexible polymer chains 
is proved Problems of influence of structural disorder on the conformational 

properties of polymer systems are widely investigated over the decades. It was shown, 
that disorder in the form of uncorrelated point-like defects of weak concentration does 
not inflict the universal properties of polymers in good solution; in particular the size 
exponent v in (HD does not change its value. This statement has been proved by 
Harris [28] and confirmed later both by renormalization group results ra and numerical 
studies [30II32] . Only when concentration of obstacles reaches a percolation threshold, 
the situation changes drastically |33ll38] : the scaling law (II|) holds with exponent 
_ Q 55 Q(^ 5 j in d = 3 [37]. Similar situation can be found in the problem of directed 
walks: some of the researches state that the value of transversal exponent changes 
to 2/3 [3^Ho] while some other appeals that it remains unchangeably 1/2 [TT] . 

Similarly as presence of external field can lead to spatial anisotropy and existence 
of two characteristic length scales of polymer macromolecule (jlj), presence of extended 
structural defects of parallel orientation (columnar defects) also causes the analogous 
effect m In this case, one should again distinguish between the components of the 
length scale in directions parallel (i?^) and perpendicular {R\) to extended defects. 
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which scale with nontrivial exponents z/y and respectively. In Ref. [12] it was shown 
that the inequality < u < z/y should hold, where z/ is a size exponent of polymer 
in pure environment ([T|). The conformational properties of flexible polymer chains in 
presence of extended defects in form of lines of parallel orientation have been studied 
in Refs. [13|IH] and numerical estimates for exponents z/_l, z/y obtained. On the other 
hand, it was shown, that the presence of extended columnar defects does not alter the 
universal properties of directed polymers [15] . 

In present work, we study the conformational properties of partially directed and 
semiflexible polymers in anisotropic environment with extended columnar defects. The 
layout of the paper is as follows. In the next Section, we introduce the lattice models 
of partially directed and biased self-avoiding walks. In Section 3, we shortly described 
the Pruned-Enriched Rosenbluth method, which is applied for analysis of properties of 
the models in Section 4. We end up with conclusions in Section 5. 

2. The Model 

To study the universal conformational properties of semiflexible and directed polymers in 
anisotropic environment, we start with the simple model of self-avoiding random walks 
on a regular lattice. We are working within the frames of Rosenbluth growing chain 
algorithm [46]: SAW trajectory grows step by step until it reaches the total length N. 
For each nth step {n < N) of growing trajectory, there are Zn available nearest neighbor 
sites to make n -|- 1th step, which variates from 0 to 2(i — 1 due to requirement that the 
SAW cannot cross itself and visit the same site more than once. If = 0 at n < A^, 
the process is stopped and a new chain is growing, until the desired total length N 
is reached. Thus, n -|- 1th step is performed towards any of Zn directions with equal 
probabilities pn = 1/^n- To generalize this model for the case of semiflexible polymers, 
we take into account that steps in directions, which lead to bending of the trajectory 
(see Fig. [T]left) are made with smaller probabilities pn • (1 — k) with stiffness parameter 
0 < A: < 1 [6]. This corresponds to biased self-avoiding walks (BSAW). With /c = 0 
one restores the simple SAW (flexible polymer chain), whereas k = 1 corresponds to the 
limit of completely stiff, rod-like polymer structures. To study the properties of partially 
directed polymers under the acting of stretching field, we take into account that the 
step with negative projection in direction x (see Fig. [1] right) could be performed with 
smaller probability p„ • (1 — p) with stretching parameter 0 < p < 1. This corresponds to 
partially directed self-avoiding walks (PDSAW). By changing p from 0 (limit of flexible 
chain - SAW) to 1 (directed polymer chain - DSAW) we can study the properties of 
conformational transitions between these states. 

We are interested to analyze the properties of these models in presence of quenched 
disorder in the form of parallel lines, stretching out along direction and randomly 
distributed in xp-plane of a lattice. Thus, the “projections” of these lines onto the xy 
plane of the lattice form a set of point-like uncorrelated defects. To realize this type 
of disorder, we start with a:p-plane of a pure lattice. Each site on this plane is then 
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Figure 1. Left: Schematic presentation of semiflexible polymer within the framework 
of BSAW model. The “gauche” steps leading to a bending of the trajectory are allowed 
with probability 1/3 ■ (1 — k). Right: Schematic presentation of partially directed 
polymer within the framework of PDSAW model. The step with negative projection 
in direction x is allowed with probability 1/3 ■ (1 —p). 


assigned to be occupied by point-like defect (and thus not allowed for SAW trajectory) 
with given probability c and empty otherwise. Then, we built the lines in direction z 
through each site containing point-like defects (see Fig. 2). We consider concentrations 
of disorder up to 40 percents due to the fact that on d = 2-dimensional lattice there is 
a percolation threshold closely to this concentration, and disorder with c > 0.4 leads to 
trivial results |44] . 

3. The Method 

The Pruned-Enrichment Rosenbluth Method (PERM) |17] is used to explore the 
peculiarities of the model. The PERM is based on original algorithm of growing 
chain proposed by Rosenbluth and Rosenbluth [46] with additional population control 
parameters. On every n-th {n < N) step of growing SAW trajectory a weight Wn is 



Figure 2. Schematic presentation of directed polymer chain in anisotropic 
environment. The longitudinal direction of polymers (x) and direction of alignment of 
extended defects (z) are perpendicular to each other. 
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assigned to a chain: 


Wn = llmi (5) 

i=l 

where mi is a nnmber of available nearest neighbor sites. Note that in calcnlating 
mi one should take into account, that in the case of BSAW or PDSAW the steps in 
directions, which lead to bending of the trajectory or steps with negative projection on 
axis X, respectively, should be properly taken with smaller probabilities (see previous 
section), which leads in general to decreasing of as comparing with Zi. Additionally, 
it is important to note that the growing trajectory cannot visit lattice sites occupied by 
defects (in our case, these sites form a set of parallel lines that are randomly distributed 
in xy plane and penetrate through the whole system). 

When the chain of total length N is constructed, the new one starts from the same 
starting point, until the desired number of chain conhgurations are obtained. Note 
that we are dealing with the case of quenched disorder, performing the conhgurational 
average of any observable over an ensemble of different realizations of disordered lattices 
according to: 


((■■■)) 



M 


M 




( 6 ) 


here is a weight of fc-th conhguration of A^-step trajectory and M is number of 
conhgurations. Though in general the critical behavior of systems with quenched and 
annealed disorder could be quite different, an equivalence of the inhuence of quenched 
and annealed disorder on scaling properties of long hexible polymer chains is proved |27j . 

While the chain grows by adding monomers, its weight huctuates. PERM 
suppresses these huctuations by pruning conhgurations with too small weights, and 
by enriching the sample with copies of high-weight conhgurations [T7jH8]. Pruning and 
enrichment are performed by choosing thresholds and hP^, which are continuously 
updated as the simulation progresses. We take [48]: = C(Z„/Zi)(c„/ci)^ and 

= 0.21T^, where c„ is the number of created chains of length n, and the parameter 
C controls the pruning-enrichment statistics; it is chosen in such way that on average 
10 chains of total length N are generated per each tour. 


4. Results and discussions 

In the present work, we study conformational properties of partially directed and 
semihexible polymers in anisotropic environment with quenched disorder in form of 
extended columnar defects. It is important to note that in our problem we have to 
perform two types of averaging of any observable: the hrst average is performed over 
all BSAW or PDSAW conhgurations on a lattice with hxed conhguration of impurities 
according to ([6|); the second average is carried out over diherent realizations (replicas) 
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Figure 3. Amplitude ratios as functions of concentration of disorder. 


of disorder: 

_ 1 ^ 

= p) 

i=l 

where M is a number of replicas. We use cubic lattices of the size L = 300 and performed 
averages over 400 replicas with a set of up to 10^ chains on each of them. 

4 . 1 . Partially directed polymers 

We analyze the properties of partially directed polymers in anisotropic environment 
within the frame of PDSAW model under varying both the stretching parameter p and 
concentration of columnar defects c. 

Let us recall, that in our model the polymers are considered to be directed along 
the x-axis, whereas the extended columnar defects are oriented along z-axis (Fig. |2]), 
so that the stretching held acts in direction perpendicular to the lines of defects. Thus, 
we have three characteristic directions in the systems, and corresponding components 
of the end-to-end distance of PDSAW are expected to behave differently: 

J^={yN- yof ~ 

W) = i^N- ZoY ~ ( 8 ) 

with non-trivial size exponents Uy, In case when no disorder is present in the 
system, we restore i^x = ^\\) (cf. (jl])). 

We start with calculation of the size exponents for completely directed polymers 
(p = 1). The data obtained lead us to conclusion that, in accordance with Ref. [45] . 
inhuence of columnar disorder is rather trivial. The values of size exponents flHl) 
governing the scaling of components of {PYiY presence of disorder are not modihed as 
compared with those of DSAW on pure lattices ((Rpure))- Thus, the anisotropy caused by 
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Figure 4. (a): Size exponents of partially directed polymers as functions of stretching 
parameter p at fixed concentration of linear defects c = 0.2. (b): Size exponent as 
function of stretching parameter p and concentration of linear defects c. 


stretching field plays the dominant role comparing with an effect of anisotropy cansed 
by disorder. At this point it is interesting to consider the amplitnde ratios 

\-^ a:,pure/ 


7 m - 7 d 2 -which are plotted on Fig. [3] as functions of concentration of disorder c. 

(-Ry.pure)’ i^z,pure) ^ ^ ^ 

These ratios may help us to catch the effects, caused by presence of extended impurities, 


even when the universality class remains unchanged. One notes that presence of parallel 
columnar defects leads to slightly shrinking of polymer in directions x and y (decreasing 
of corresponding amplitude ratios) and correspondingly to slight expansion in direction 


2;. 


To find the numerical estimates for the critical exponents in (|H]), we analyzed the 
behavior of corresponding components of the end-to-end distance as functions of the 
length of growing PDSAW trajectory under varying both the stretching parameter p 
and concentration of defects c and applied the least-square fitting to the data obtained. 
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The results obtained are presented on Fig. 0^. Let us recall, that the critical exponents 
of simple SAW are universal in the sense that they are dependent on the only one “global 
parameter” - the space dimension d (see e.g. Eq. ([T])). In our case of PDSAW, we note 
the non-universal behavior of exponents u as functions of the stretching parameter 
p, which now plays the role of another global parameter. Note, that similar non- 
universality is obtained in renormalization group studies of BSAW model in Ref. 09] , 
where, in particular, the dependence of critical exponent u on the stiffness parameter 
(bending energy) is obtained. If the stretching field is absent {p = 0), we restore 
the problem of simple SAW in presence of linear defects studied in jUj. In this case, 
as expected, one hnds stretching of trajectory in direction, parallel do orientation of 
defects, so that > Vx = On the other hand, the competition between the 
stretching held and space anisotropy caused by presence of defects leads to existing 
of three characteristic length scales in the system. With increasing the parameter p 
at hxed concentration of disorder, the stretching of polymer in direction x increases 
{i/z > Vx > I'y) and hnally overcomes the effect caused by disorder {vx > 
for large values of p). Thus, at each hxed concentration of disorder, there should 
exist some critical value of stretching parameter p, when the anisotropy caused by 
presence of extended defects is exactly “counterbalanced” by the stretching held, acting 
in perpendicular direction, so that Ux = Vz = ^xz- This may be treated as a transition 
point between the state where the role of disorder is dominant and the polymer is mainly 
oriented along z axis, and the state when the stretching held starts to dominate and 
causes the polymer chain reorientation mainly along x axis. Analyzing our data under 
varying both the stretching parameter p and concentration of columnar defects c, we 
found the pairs of c and p values corresponding to counterbalance between two types of 
anisotropy, and estimated the values of exponents presented on Fig. |Dd. 

4-2. Semiflexible polymers 

We analyze the properties of semihexible polymers in anisotropic environment within the 
frames of BSAW model under varying both the stihness parameter k and concentration 
of columnar defects c. Numerical simulation are performed for the chain lengths up to 
100 monomers. 

We start with evaluating the number of bends in a polymer chain of hxed length 
N, corresponding to a number of times when BSAW trajectory changes its direction, 
as function of parameters k and c (Fig. |5^). As expected, nf, decreases gradually 
with increasing the stihness of a chain. At small values of fc, corresponding to regime 
of hexible polymer chain, presence of columnar defects does not alter the bending of 
trajectory, whereas for considerably stih chains (with k close to 1) the number of bends 
slightly decreases with increasing the concentration of disorder. This ehect can be 
visualized more clearly by analyzing the average number of steps between two successive 
bends in BSAW trajectory (average length of a straight segment Igs, which can be 
related to the persistence length of semihexible polymer), which increases with disorder 
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concentration (see fig. Eb)- For fixed chain length we can analyze the behavior of Igs as 
function of 1 — k (see £g. ED. These numerical data can be approximated according to: 

hs = a(l - kf, (9) 

where the values of constants a and b depend on the concentration of columnar disorder 
(see Table [ID. 

Note, that one can introduce a bending energy e associated to each “gauche” step, 
so that the statistical weight e“^ corresponds to each bend of the BSAW trajectory. 
With e: = 0 we restore the flexible polymer chain {k = 0), whereas with increasing £ 
the bends become energetically unfavorable, so that e ^ oo corresponds to limit of stiff 
rod-like polymers {k = 1). One thus has [6]: 

1 — k = e~^. (10) 

The total energy of Wmonomer semiflexible polymer chain is thus given by = Ub-e. 
Statistical fluctuations of the energy indicates transitions or crossovers between 




Figure 5. (a): Number of bends in BSAW trajectory of total length N = 100 at 
different concentrations of disorder c as function of stiffness parameter k. (b): Average 
length of straight segment as function of chain length for different concentrations 
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Figure 6. Average length of straight segment as a function of stretching parameter 
for different concentrations 

different conformational states. In the case of semiflexible polymers, this corresponds 
to the transition between regimes of flexible and stiff chains. We considered generalized 
“specific heat” per one monomer defined via energy flnctnations as: 

(IM)-W^), ( 11 ) 

and analyzed the peak strnctnre oi Cn as fnnction of e at varions concentration c of 
linear defects (Fig. [7]). In the case of pnre environment, we hnd for the critical valne 
of bending energy £ 1.47 (corresponding to k = 0.77). This valne is shifted with 

increasing the concentration of disorder (at c = 0.2 we observe £ ~ 1.24 {k ~ 0.71). 
Thns, as expected, the presence of extended impnrities makes the polymer chain stiffer 
and canses the conformational transition into the rod-like state to take place at smaller 
valnes of bending energy e: comparing with the case of pnre environment. 


c 

a 

b 

0 

1.234(6) 

-0.382(5) 

0.1 

1.229(6) 

-0.385(6) 

0.2 

1.224(6) 

-0.394(6) 

0.3 

1.217(7) 

-0.411(6) 

0.4 

1.205(8) 

-0.444(8) 


Table 1 . Values of parameters a and b in Eq. at different concentrations of disorder 
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Figure 7. Specific heat Cn as function of bending energy at different concentrations 
of disorder. 

5. Conclusions 

The aim of the present paper was to analyze the conformational transformations in 
polymers in solutions in presence of colnmnar obstacles of parallel orientation, causing 
the spatial anisotropy. We consider two classes of polymers: the so-called partially 
directed polymers with preferred orientation along direction of the external stretching 
held and semihexible polymers. 

We are working within the frames of model of self-avoiding random walks (SAWs) 
on a regular lattice. To study the properties of partially directed polymers, we introduce 
the stretching parameter 0 < p < 1, and by changing p from 0 (limit of hexible 
chain - SAW) to 1 (directed polymer chain - DSAW) we analyze the conformational 
transitions between these states. This corresponds to partially directed self-avoiding 
walks (PDSAW). To generalize the SAW model for the case of semihexible polymers, we 
introduce the stihness parameter 0 < fc < 1: with k = 0 one restores the simple SAW 
(hexible polymer chain), whereas k = 1 corresponds to the limit of completely stih, 
rod-like polymer structures. This corresponds to biased self-avoiding walks (BSAW). 

Our numerical analysis of PDSAW model in an anisotropic environment with 
extended columnar defects reveals a new universality class with non-trivial values of 
critical exponents. The competition between the stretching held and space anisotropy 
caused by presence of defects leads to existing of three characteristic length scales in 
the system. At each hxed concentration of disorder c, we found the critical value of 
stretching parameter p, when the anisotropy caused by presence of extended defects 
is exactly “counterbalanced” by the stretching held, acting in perpendicular direction. 
This may be treated as a transition point from the conformation when the polymer is 
mainly oriented along the position of extended linear defects to the regime when the 
stretching held starts to dominate and make the polymer chain reorientation mainly 
along the perpendicular direction. 
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Numerical simulations of semiflexible polymers in anisotropic environment also 
reveal difference in comparison with pure environment. As expected, presence of 
extended columnar defects of parallel orientation leads to an increase of the stiffness of a 
chain. In particular, the average number of steps between two successive bends in BSAW 
trajectory (average length of a straight segment Iss, which is related to the persistence 
length of semiflexible polymer) increases with disorder concentration. Analysis of the 
peak structure of specihc heat shows the shift of bending energy value with increasing 
the concentration of disorder and thus, the presence of extended impurities causes the 
conformational transition into the rod-like state at smaller values of e comparing with 
the case of pure environment. 
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